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Abstract 

We show how the basic operations of quantum computing can be expressed 
and manipulated in a clear and concise fashion using a multiparticle version 
of geometric (aka Clifford) algebra. This algebra encompasses the product 
operator formalism of NMR spectroscopy, and hence its notation leads directly 
to implementations of these operations via NMR pulse sequences. 
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Geometric (aka Clifford) algebra is a generalization of vector algebra to arbitrary dimen- 
sions and signatures, which provides a concise and geometrically transparent notation for 
describing a wide range of physical phenomena (for introductions and examples, see [|1]-|6|]). 
Multiparticle geometric algebra is a recent further generalization that enables one to deal 
with interacting two state quantum systems J?H§. Since most models of a quantum com- 
puter are based on such systems P-pH, it is of interest to formulate the basic operations of 
quantum computing in these terms, with the goal of gaining deeper insight into them. Given 
the recently discovered methods of emulating a quantum computer via NMR spectroscopy 
]T2]^H[ , it is further of interest to note that one of the main tools that NMR spectroscopists 
have developed to aid them in understanding their experiments, known as the product op- 
erator formalism H15-18H, is a subalgebra of a multiparticle geometric algebra. Geometric 



algebra also encompasses the quaternion methods often used by NMR spectroscopists to 
calculate the effects of composite pulses on product operators [p^ -plH- Thus, as illustrated 
in this letter, it is generally straightforward to translate a quantum logic operation expressed 
in the multiparticle geometric algebra into an NMR implementation thereof. 

We shall begin with a physically motivated definition of the geometric algebra <?(3) 
of a single particle. This algebra is isomorphic to the Pauli matrix algebra, where the 
isomorphism is obtained by identifying the Pauli matrices ^_i,a 2 ,^.s with an orthonormal 
basis (Ti, cr 2 , er 3 of a Euclidean vector space. Note that the product of all three Pauli matrices 
<Zli<Zl2<2l3 is j us t the imaginary unit "z" times the identity matrix, which we denote by 
l = il. This enables us to further identify the corresponding element t = <j\ <t 2 <t 3 of the 
geometric algebra £7(3), which is commonly called the unit pseudo-scalar, with the imaginary 
unit itself. In this way £(3) becomes an algebra over the real numbers, even though the 
Pauli matrix algebra is usually regarded as a complex algebra. 

The (nonrelativistic) iV-particle geometric algebra Q N (3) consists of a direct product of 
N copies of (?(3) 

g N (3) = 0(3) ® • • • <g> 0(3) (N factors) (1) 

Throughout this letter, we shall assume that 1 < i, j, k < 3, 1 < £, m, n < N, and that 
a, /3, . . . are real scalars. Then the multiparticle geometric algebra can also be defined by 
the following five straightforward rules: 

{(Ti} is a basis for a real vector space; (2) 

= ; (3) 

o\{atr? + = aoitr? + Pofrl and (4) 

{aoi + p<r?)<rl = aa\a\ + /3« ; 



25ij (the Kronecker delta); (5) 



for all t^m: tr\crf - a™(r\ = . (6) 

The NMR product operator formalism relies upon a large repertoire of relatively complicated 
rules for predicting the evolution of the spins, which vary with the interaction (e.g. strong 



This definition can be derived from a relativistic multiparticle theory based on the geometric 
algebra of N copies of space-time 
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or weak coupling), the total angular momentum (for spin > ~), and so on. With a little 
practice, however, all the rules of the product operator formalism can be readily derived as 
they are needed from the five rules given in Eqs. (|2|-§]). 

By our definition of G(3) above, a faithful matrix representation of this (2 3 )-dimensional 
real algebra is obtained by mapping the basis vectors to the Pauli matrices: cr fc — > g_ k . A 
matrix representation of Q N (3) is obtained simply by taking iV-fold Kronecker products of 
these matrices. Since the representation by Pauli matrices is a complex algebra, however, 
its real dimension is 2 ■ 4^ = 2 2N+1 , which for iV > 1 is less than the dimension 2 3N of 
Q N (3). It follows that the direct product representation is no longer faithful. The extra 
degrees of freedom are due to the fact that, unlike the matrix representation, there is a 
different complex unit { = a{ <t\ a\ for every particle. Since the matrix algebra is believed 
to contain all possible quantum states, operators and propagators, these degrees of freedom 
have no known physical relevance. They can easily be removed by multiplying through by 
a primitive idempotent of the form 

C = 1(1 - t h 2 ) |(1 - Lh 3 ) ■ " • |(1 - Lh N ) . (7) 

This correlator commutes with everything in the product algebra, and projects it onto an 
ideal of the correct dimension. The projection can be interpreted physically as locking the 
phases of the various particles together. For this reason, we shall not specify which factor 
our imaginary unit i comes from, since all choices are rendered equal by the correlator - 
whose presence in all our expressions will also not, in the interests of brevity, be written out 
explicitly. Further discussion of these issues may be found in [0,|J . 
Another class of primitive idempotents that we shall need are 

E£ = |(1 ±<) (8) 

and products thereof from different factors. These are easily shown to have the following 
properties: 

(E™) 2 = E™ , E™E™ = , E™ + E m = 1 , <r™E™ = ± E£ . (9) 

A product of N such idempotents, one for each particle, has the (2^) x (2^) matrix repre- 
sentation 

E £1 ®---®E ejv <-> \e 1 ---e N )(e 1 ---e N \ , (10) 

where e m = ±1 and E ± = |(1_ ± <r 3 ). To illustrate the utility of these idempotents, we 
will begin by using them to parametrize an arbitrary entity from the multiparticle geometric 
algebra £7^(3) in terms of "rotors" from its even subalgebra 0. We shall do this for only 
two particles, from which the general case should be clear. 

Given an arbitrary multivector Mg^(3)®^(3), we may write 

M = M(E^ + E^)(E 2 h + E 2 _) = ME}E 2 + + -- + ME^. (11) 
Let (M) + = |(M+M) denote the projection of M onto the even subalgebra (£(3)<g>£(3)) + = 
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£ + (3) <g>£ + (3), where M is the grade involution of M.Q Then, using the fact that E™E™ = 
and E™ = E™, each term above may be written as 



MEJX = (ME^^E 2 

= (ME £ \E 2 2 + ME^Ei J(E^E£) (12) 
= 2(ME ei E 2 2 > + (E ei E 2 2 ) 

(ei,e 2 G {±1}). It follows that 

M = (V> ++ E^ + V+_E^E 2 _ + ^_ + E\E 2 + + V— E^E 2 ) , (13) 

where each rotor t/> e , e „ is in £ + (3) ® £ + (3), and each ^E^E 2 , can be viewed two- 
particle spinor J7IJ9]. A similar expansion holds for any number of particles and, when 
translated into matrices via Eq. fllOD above, corresponds to an expansion of the matrix 
into a sum of matrices each of which contains one column from the original matrix, and is 
otherwise zero. 

If we multiply M by its reverse M, we obtain 

MM = (V ++ E^E^ + --- + V__E^E 2 _)(E^E^ ++ + --- + E^E 2 _^__) (14) 
= ^ ++ E\E 2 + ^ ++ + V+-E^E 2 + V-+ELE 2 + V>_ + + ^__E l _E 2 _^__ ■ 

Since the reverse corresponds to the Hermitian conjugate of matrices, this expression repre- 
sents a general Hermitian matrix. A rank one Hermitian matrix is obtained from any single 
term of the expression, e.g. V ++E 1 l _E 2 h i/> ++ . Provided that = 1, the scalar part 

of this expression is ~, which corresponds to a matrix trace of one. Thus this expression 
represents the density matrix of a pure state, while a general density matrix is a convex 
combination of terms as in Eq. (|T4T ) above. 

A simple rotation by an angle a about an axis a = SjajCTj (||a|| = 1) in a one-particle 
space is given by the exponential of its generator, e.g. exp(— ctta/2), just as is commonly 
done with Pauli matrices. Due to the commutivity of vectors cr™, cr™ from different particle 
spaces, one cannot express arbitrary rotations between particle spaces in the algebra. An 
exception to this rule are the particle interchange operators, which in the case of two particles 
is given by 

H 1 ' 2 = |(1 + <t\o\ + <t\<t\ + <r\*l) . (15) 

An important variation on this is given by the scalar coupling Hamiltonian in NMR spec- 
troscopy, J 1 ' 2 = 7rJ 1,2 (2II ' 2 — l)/2, where J 1,2 is a coupling constant in Hertz [17]. Since it 



is easily shown that II ' 2 squares to one, the corresponding time- dependent propagator is 
exp(tJ 1 - 2 t) = e-^ J1,2 * /2 (cos(7rJ 1 - 2 t) +tn 1 ' 2 sin(7rJ 1 ' 2 t)) . (16) 



2 For definitions of the grade involution, the reverse, and other common geometric algebra terms, 
see §. The identity (G(3) ® <?(3))+ = G + (3) ® G + (3) holds because for all 1 < i,j < 3 and 
1 < m, n < N, er™<x™C = — i m cr™i n cr™C where C is the correlator defined in Eq. (|7|). 
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In the product operator formalism, one generally assumes weak coupling, i.e. \uj 1 — uo 2 \ ^> 
7r J 1 ' 2 , where u 1 , u 2 are the resonance frequencies of the spins. This enables the usual two-spin 
Hamiltonian of NMR spectroscopy, i.e. H 1 ' 2 = J 1 ' 2 + K 1 ' 2 where K 1,2 = ^{u l cr\ + uj 2 cr\) 
denotes the Zeeman Hamiltonian, to be replaced by its first-order approximation 



H^ 2 ak ss K 1,2 + \^a\a\ . (17) 

Since this approximation is diagonal in the usual cr\ + a\ basis, its propagator can be 
written down in closed form (see e.g. |H)|-|T7|). Geometric algebra is however not limited to 
weak coupling, although strong coupling of course complicates the propagators substantially. 
This will be demonstrated shortly, when we show how to diagonalize a general two-spin 
Hamiltonian. 

We next show how to implement a rotation in one particle space conditional on the state 
of another it is correlated with. The simplest example is the controlled-NOT (or XOR) 
quantum gate, which is generated by a transition Hamiltonian of the form 

H t ? == ^E 2 _ = \a\{\-al). (18) 

Since {cr\) 2 = 1 and cr\ commutes with the idempotent E 2 ,, the corresponding propagator 
is given by 

nf{a) = e" taH - 2 / 2 = e^/ 2 E 2 +E 2 ,. (19) 

This rotates the first spin about cr\ in those states in which the second spin is along the 
— <Xo axis. The controlled-NOT n\ 12 ee '(it) is obtained when a = ir, in which case the 
density matrix of the basis states E^E 2 2 are transformed as follows: 



R!' 2 (E e \Ei)R!' 2 = K„K = { e^e! It : l\ < 20 > 
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In order to obtain a controlled-NOT that preserves the phases of the basis states (as is 
usually assumed), one need only multiply R^' 2 by the conditional phase shift exp(*/7rE 2 /2) = 
lEi 2 _ + E^. Similarly, the propagator for the controlled-controlled-NOT (or Toffoli gate) can 
be written as 

e «r(l-^)E» E»_/2 = a l E 2_ E 3 + ^ _ g2 ) ^ (21) 

while the Fredkin gate is given by the conditional particle interchange 

e «r(l-n^)E3 /2 = n l,2 E 3 + E 3 (22) 



In , implementations of the controlled-NOT and Toffoli gates by NMR pulse sequences 
were given, which can also be easily validated by geometric algebra methods. In the case of 
the controlled-NOT, the implementation consists of a product of three propagators, which 
may be simplified as 
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S^' 2 = e ~ CTCT i/ 4 e - OT0 3 £r l/ 4 e - wro 2/ 4 

= e~ t ^ CT i/ 4 (e~ twCT 3 CT ^ 4 e~ t ^2/ 4 e ^ CT 3 cr 3/ 4 )e _i ' 7rcr 3 <T l/ 4 (23) 

_ e -t.7r -l/4 e i,7r<Tj C r|/4 e -i,7r<Tlo-2/4^ E ^ + g2 ^ 

= (e-^i^El + E 2 f )(e- tmT 3/ 4 E 2 h + e ^/ 4 E 2 _) 



where we have used Eq. @ to get the relation exp(— t7rcr3cr|/4)E 2 t = exp(=Ft7rcr3/4)E 2 t . To 

get a controlled-NOT that preserves phases, we may left-multiply S^' 2 by exp(— nr(l + cr\ — 
<r§)/4) as in H, or right-multiply it by exp(— tvrcrgE 2 /4) exp(t7rcr^E 2 h /4) = exp(t7T(T3cr|/4) 
along with exp(t7rE 2 _/2) as above. The validation of the Toffoli pulse sequence is similar, 
though considerably more complex. 

An interesting generalization of the two particle conditional rotations is obtained by 
taking the exponential of all four transitions, i.e. 

exp(-/,(a< + 0-}E 2 + c^o-jE^ + /3+<r 2 E|_ + /3_<r 2 Ei)/2) 
= exp(-t(( a+ + a-)tr{/4 + (J3+ - /3_)^o- 2 /4)) (24) 

exp(-t((a+ - a-)a\ai/A + (0+ + P-)o-\/A)) 
= exp(— tX) exp(— lY) = exp(— LY) exp(— tX) , 

where X and Y commute. In addition, since 

X 2 = (( a+ + a_) 2 + (/3 + -/3_) 2 )/16 (25) 
Y 2 = (K-«_) 2 + (/3 + +/5_) 2 )/16, 

their propagators can be written in closed form as 

exp(-tX) = cos(v / X 2 ) - tsin(v / X 2 )X/v / X 2 " (26) 
exp(-tY) = cos(v / Y 2 ) - tsin(v / Y2)Y/v / Y 2 " . 

For example, if a + = (3 + = and «_ = /?_ = 7ta/2, we obtain (0-3 + erf; — a\o\ — cr\a\)/2, 
which is the particle interchange operator II ,2 up to a conditional phase of exp(t7rE 1 |_E 2 F ). 
If, on the other hand, a + = (3- = 7ta/2, we get (c^ — <r| — ct\(t\ + <j\(t'\) /2, which is converted 
by exp(t7rE ] | _E 2 ) into |(1 — cr\cr1 + cr\o\ — cr\cr!£). The terms 1 — cr\cr1 are annihilated by 
E^E^, while cr\cr\ — a\cr\ induces direct transitions between the E+E 2 ^ and E^E 2 states: 



12 12 
(T 1 cr 1 — er 2 er 2 



lE^E|(o-}o- 2 - ct\(tI) (27) 
= E*JS*(«ri«r? - via*? = E^E|(2 + 2<7 3 <7 3 ') = 4E^E 2 F . 

We note that these transformations can be implemented in NMR via what might be called 
compound pulses, i.e. multiple simultaneous pulses each selective for a single transition 

Conditional rotations further enable us to diagonalize a general two-spin NMR Hamil- 
tonian H 1 ' 2 as above. If we define the correlated idempotents E± 2 = |(1 ± erger 2 ), this may 
be rewritten as 
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H 1 ' 2 = y<rl + iaj 2 <r 3 2 + tt J 1 ' 2 + <rJ<r?E^ a ) . (28) 

We can split H 1 ' 2 into two commuting parts as follows: 

H 1 2 = H 1 ' 2 E^ 2 + H 1 - 2 E^ 2 = H^ 2 + H^ 2 . (29) 

Since E^Ei' 2 = 0, the H 1 ^ 2 part is already diagonal. The off-diagonal part, on the other 
hand, can be written as 

Hi: 2 - ivrJ^Ei: 2 = (u-ffl + tvJ^ctIctDE 1 : 2 , (30) 

where u± = uj 1 ±uj 2 . The 0-3 and cr\cr\ terms are rotated into one another by ex~p((f)Lcr\cr\) , 
and hence this can be written in polar form as 

H^ 2 - ±vr J^Eh 2 = e - 3 exp(^ -> 2 )E^ , (31) 

where G = \J iyj^) 2 + (vrJ 1 ' 2 ) 2 and = arctan(7r J 1,2 /u;_). The right-hand side of this 
equation is clearly diagonalized by exp(— (f)L(r\cr\/2) and commutes with El_' 2 , from which 
it follows that H 1 ' 2 is diagonalized by the conditional rotation 

T(0) = e -^^B-' a /2 = e -^i/2 E h2 + E l,2 (32) 

The result is a weak coupling Hamiltonian whose frequencies have been shifted by ±0, 

f (0)H 1 - 2 T(0) = + 6)<x 3 + - 6)<x 2 + IvrJ 1 ^^ 2 , (33) 

while the corresponding transition moment shows that the peak intensities have been altered 
by 

(KTiM'l + *?W)Ey [E 2 2 f (0)(<r{ + <r 2 )T(0)E 2 j) (34) 
= (1 + e! sin(0))5 eii£ /5 £2i _ e / + (1 - e 2 sin(0))5 ei _ £ ;<5 e2 , ei . 

The Hadamard transform plays an essential role in many quantum algorithms, but its 
simple geometric interpretation is seldom pointed out. Consider a rotation of the m-th 
particle by an angle a about the w m = (cr™ + a™)/\/2 axis: 

W m (a) = exp(-taw m /2) = cos(a/2) - iw m sin(a/2) (35) 

The one-particle Hadamard transform W m = —iw m is obtained when a = n, and the 
iV-particle Hadamard transform W^r is simply the commutative product of the Hadamard 
transforms of the individual particles. The quantum Fourier transform (QFT) |23[ can 
likewise be written, and simplified, using the multiparticle geometric algebra. As shown by 
Don Coppersmith (unpublished manuscript), the iV-particle QFT Qa? can be written as a 
recursive product of one-particle Hadamard transforms and two-particle conditional phase 
shifts, 

V £ ' m = exp(-tu em E e + E™) (uj em = ir2 e - m for I < m) , (36) 
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namely 

Q N = U 1 --.^- 1 ^ = (W 1 V 1 ' 2 ---V 1 ' 7V )---(W JV - 1 V^ 1 ' iV )(W Ar ) . (37) 

Each factor U m above can be rearranged so that all the Hadamaard transformations 
come first: 



U m = (-tw m )(e-^ m+lE + E + +1 ---e- l ^ E + E ?)(tw m )(-tw m ) 

_ e - t (- t w m )E™( f .w m )(o; m , m+1 E™ +1 +-+a; miiv E^)^_ tw m^ /gg-j 
_ e - t /2 e - tC r™(a; m , m+ iE™+ 1 +-+^ m ,ivEy) ^_ tw m-j 



Since the factors U n with n > m do not involve the m-th particle space, we may therefore 
write the QFT with all the — iw m together on one side, i.e. 

Q N = e -W 2 e - tCT iK2E2 + ... + o, liJV EN) ... e - lCT f- 1 (^_ 1 , JV E^) Wjv _ (gg) 

In NMR spectroscopy, each of the conditional rotations in this expression can in principle be 
implemented with a single compound pulse. The Hadamard transform Wjv of all N spins, 
on the other hand, can be implemented with three "hard" (nonselective) pulses, namely 
exp(t7r<T2/8) exp(t7rcr 1 /2) exp(— 17:02/8), which takes only a very small fraction of the time 
required for a "soft" pulse (or the conditional rotations). The complexity of this NMR 
implementation of the QFT is only O(N), and since the time required for a hard pulse is 
essentially independent of the number of spins involved, the implementation requires only 
half the number of steps in the parallel implementation proposed by Coppersmith. 

In closing, we note one final advantage of multiparticle geometric algebra, which is the 
astonishing ease with which one can compute partial traces. In terms of idempotents, the 
partial trace of a density matrix p over any one particle m can be written as: 

Tr m (p) = E™pE™ + E^pE™ + <rf (E™pE™ + E™p E™)< . (40) 

Consider, for example, the three particle GHZ state, Pghz = 

(1 + o\o\ + o\o\ + o\o\ + o\o\o\ - v\(t\<j\ - (t\ct\<t\ - o\o\a%)l& . (41) 

Using the relations Ej_cr|E| = E|E|cr| = and cr|E|cr|, = E| for k = 1,2, one finds that 

E|p GHZ E 3 ± = E%{1 + trio* + trl** + ol<rl)/8 (42) 

and 

<t?E 3 ±Pghz E|<x? = E|(l + olo\ - o\o\ - o 2 3 ol)/8 . (43) 

The sum of the four terms in Eq. ( f4"0"D is thus the mixed state 

Tr 3 (p GH z) = (1 + <^3)/4 . (44) 

More generally, the partial trace is obtained simply by dropping all terms from the density 
matrix containing factors from the particles over which the trace is taken, and multiplying 
by two to the number of such particles. 
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In conclusion, we have shown how to formulate the most important operations of quan- 
tum computing in multiparticle geometric algebra, and illustrated some of the advantages 
that this more general theory has over the established product operator formalism used in 
NMR spectroscopy. These advantages will be further demonstrated in subsequent publica- 
tions devoted to developing the theory and methods needed for ensemble quantum computing 
by NMR spectroscopy. 
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